Experimental characterization of bosonic photon creation and annihilation operators 
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The photon creation and annihilation operators are cornerstones of the quantum description of 
the electromagnetic field. They signify the isomorphism of the optical Hilbert space to that of 
the harmonic oscillator and the bosonic nature of photons. We perform complete experimental 
characterization (quantum process tomography) of these operators. By measuring their effect on 
coherent states, we obtain their process tensor in the Fock basis, which explicitly shows the "raising" 
and "lowering" properties of these operators with respect to photon number states. This is the first 
experimental demonstration of complete tomography of non-deterministic quantum processes. 



Quantum operators of annihilation a and creation a 1 
of bosonic particles act on states with a definite number 
m of identical particles, resulting in that number being 
incremented or decremented, respectively: 



|m) = vm + 1 \m + 1) 



(1) 
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First proposed by Dirac in his 1927 Quantum Theory of 
the Emission and Absorption of Radiation [l\ . these oper- 
ators play an enormous role in many fields of physics and 
chemistry: quantum mechanics, quantum optics, quan- 
tum chemistry, quantum field theory and condensed mat- 
ter physics. Specializing to optics, they are instrumen- 
tal in quantum description of light, giving rise to many 
fundamental phenomena such as spontaneous emission, 
Lamb shift, Casimir force and lasing. Equally impor- 
tant is practical implementation of photon creation and 
annihilation, which provides us with a universal toolbox 
for manufacturing arbitrary quantum states of light, re- 
quired for quantum information processing and quantum 
communications 0, El ■ 

Implementation of a and d^ is however challenging. 
This is because these operators do not preserve the trace 
of a state's density matrix, which means they cannot oc- 
cur in the framework of deterministic Hamiltonian evo- 
lution of a physical system. Therefore bosonic creation 
and annihilation can be obtained in the laboratory only 
in an approximate, non-deterministic fashion. That is, 
the action of operators occurs with probability less than 
one, but is heralded by a classical event. 

First successful realization of the optical photon cre- 
ation operator in this manner in 2004 Q gave rise to a 
new class of states known as photon added states. In 
2006, Ourjoumtsev and colleagues applied the photon 
annihilation operation to the squeezed states generat- 
ing optical "Schrodinger cats" [5|. Neergaard- Nielsen et 
al. generalized this approach in 2010 to generating arbi- 
trary continuous- variable qubits @. In 2007, Parigi et 
al. verified non-commutativity of a and a* in application 



to the thermal state Q- The photon annihilation oper- 
ator has been used for continuous-variable entanglement 
distillation in 2010 by Takahashi and co-workers [8|. Ex- 
perimental recording of photon creation and annihilation 
events in the time domain has been reported by Gleyzes 
et al. 0. 

In order to include the photon creation and annihila- 
tion operations into the quantum technology toolbox, we 
need to develop methods of their characterization and 
performance testing. This is the purpose of the present 
work. We implement a and experimentally and ana- 
lyze them as quantum "black boxes" , or quantum pro- 
cesses. By probing them with coherent optical states 
(weak laser pulses) of different amplitudes and measuring 
the quantum fluctuations of the output electromagnetic 
field, we determine how these black boxes would affect 
any arbitrary state of light within a practically relevant 
subspace of the optical Hilbert space. As a result, for the 
first time since the photon creation and annihilation op- 
erators have been discovered, we explicitly observe their 
action on the photon number states to be in accordance 
with Eq. flTJ. 

The method we employ for the characterization of 
quantum processes relies on the optical equivalence the- 
orem. According to that theorem, the density operator 
p of an arbitrary state can be written as a linear combi- 
nation of coherent-state density operators, 



P 



J Pp(a) I a) ("I 



d a, 



(2) 



where Pp(a) is the Glauber-Sudarshan P function of state 
p. Further, since every quantum process £ (in this case, 
photon creation and annihilation) is a linear map with re- 
spect to density matrices, we can write the process output 
as 



£(p) = J Pfi£(\a) (a\)d 2 a. 



(3) 



If we know £ (|a) (a\) for every coherent state |a), we can 
determine the process output £{p) for any state p. 
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This is useful because optical states that are employed 
in quantum information processing (for example, number 
states or their superpositions) are highly nonclassical and 
cannot be generated easily. In contrast, coherent states 
are directly obtained from lasers. Our method permits 
us, by probing the "black box" with simple laser pulses, 
to learn its effect on any other state of light, however 
complicated it may be. In the past, this approach, re- 
ferred to as coherent-state quantum process tomography 
(csQPT) fiol 12 1 . has been applied to the processes of 
attenuation, phase shift [l(J and quantum optical mem- 
ory ■ A closely related method has been used for the 
quantum characterization of optical detectors 3, HI 

The result of csQPT — the data about the process 
- can be compactly written in the form of a process 
tensor. This is a rank-4 tensor £V^ n such that, for any 
process input p, the density operator of the process out- 
put in the photon number basis is given by [£{p)]- k — 
T m n £jk n P mn - The process tensor is calculated accord- 
ing to 



a) 



£T= / P mn (a)(j\£(\a)(a\)\k)d 2 a, 



(4) 



where P m n(ct) is the P function of operator |m) (n|. Com- 
putation of the process tensor is complicated by highly 
singular nature of this function; Refs 13-12, [HI elab- 
orate different ways of resolving this complication. An- 
other practical issue is associated with the infinite dimen- 
sion of the optical Hilbert space. In csQPT experiments, 
the process tensor is evaluated for a subspace H(n max ) 
spanned by number states up to a certain cut-off value, 
n max . The choice of n max is determined by the maximum 
amplitude a max of the set of coherent probe states used 
in the experiment, as well as the reconstruction method 
used. In our work, n max = 7. 

Practical realization of photon annihilation employs a 
low reflectivity beam splitter, through which the target 
state is transmitted [Fig.QJa)]. Detection of a single 
photon in the reflection channel indicates that this pho- 
ton has been removed from state In this case, the 
state emerging in the transmission channel of the beam 
splitter is approximated by d\ip). 

For photon creation, low-amplitude spontaneous para- 
metric down-conversion (SPDC) in a nonlinear optical 
crystal can be employed. The target state enters the sig- 
nal SPDC mode from the back of the crystal [Fig. [Tfb)] . 
If SPDC occurs, operators a) act simultaneously on both 
the signal and idler modes. If this event is heralded by 
the photon detector in the idler channel, the target state 



0- 



becomes a 1 

In the experiment, we employ a mode-locked 
TkSapphire laser (Coherent Mira 900), which emits 
transform- limited pulses at ~790 nm with a repetition 
rate of 76 MHz and a pulse width of ~1.6 ps. Herald- 
ing photons are registered by PerkinElmer SPCM-AQR- 
14-FC single-photon detectors. The field quadratures of 
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FIG. 1. Experimental setups for photon annihilation (a) and 
creation (b). The processes are heralded by "clicks" in single- 
photon detectors. 



output states are measured by means of high-bandwidth 
balanced homodyne detectors 



18 20]. Both the probe 



field and the local oscillator field for homodyne detection 
are obtained from the master laser. The amplitude of 
the probe field is varied using a half-wave plate and a 
polarizing beam splitter. 

In order to obtain SPDC, required for the photon cre- 
ation, the light from the master laser is frequency dou- 
bled in a single pass through a 17-mm long lithium trib- 
orate (LBO) crystal, yielding a typical ^80 mW average 
second-harmonic power after spatial filtering. This field 
is focused, with a waist of 100 ^m, into a 2-mm long 
periodically poled potassisum-titanyl phosphate crystal. 
This crystal is phase-matched for type II SPDC, with 
the signal and idler modes being spatially and spectrally 
degenerate but having orthogonal polarizations. 

We acquire a set of field quadrature data for the out- 
puts of both processes for a set of probe coherent states 
with amplitudes ranging from to about 1.7. For the 
photon annihilation process, the output states are al- 
most identical to the input states, as they only undergo 
slight attenuation propagating through the beam split- 
ter [Fig. [Ha)]. This is to be expected, because coherent 
states are eigenstates of a . Photon creation, on the other 
hand, strongly changes the nature of the state, produc- 
ing single-photon added coherent states, studied in de- 
tail theoretically by Agarwal and Tara 21|, experimen- 
tally by Zavatta et al. [3, Q • Highly nonclassical nature 



FIG. 2. The Wigner functions of coherent states of different 
amplitudes subjected to the photon creation operator. 



of these states is evidenced by their Wigner functions 
(Fig. [2]), which are explicitly non-Gaussian and negative 
in some regions of the phase space. 

The processes we study are non-deterministic, and 
their probability of occurrence depends on the input 
state. Accounting for this dependence is crucial for the 
correct reconstruction. In csQPT, this is done by renor- 
malizing the process output for the probe states so that 
Tr[£ (|a) (a\)] in Eqs. ((3]) and (j4|) is proportional to the 
probability of the heralding event [llj . To illustrate 
the significance of this step, it is instructive to apply 
Eq. ([3]) to the photon annihilation operator, such that 
£(\a) (a\) = \a\ 2 \a) (a\. If the coefficient \a\ 2 , responsi- 
ble for the non-deterministic nature of a, is neglected, we 
would obtain the identity process. We see, remarkably, 
that the "lowering" feature of a arises in csQPT entirely 
due to the variation of the event probability as a func- 
tion of the probe amplitude, rather than transformation 
of the probe state itself. 

The information on the heralding event probability is 
obtained by keeping track of the photon count rates for 
various input states. Theoretically, we expect these rates 
to behave as 

pr a (a) oc (a| a a |a) = a 2 ; pr a t (a) oc (a| aa^ \a) = 1+a 2 . 

(5) 

The experimentally observed dependencies are consistent 
with these expectations as displayed in Fig. [3l 

We use the iterative algorithm [l2| to reconstruct the 
process tensors directly from the acquired field quadra- 
ture data. The algorithm makes use of the Jami- 
olkowski isomorphism between quantum processes ap- 
plied to Hilbert space % and positive semidefinite oper- 
ators over the Hilbert space Hfgi'H [22[ . In this way, the 
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FIG. 3. Event count rate for the photon annihilation (a) and 
creation (b) operators as a function of the input coherent 
state amplitudes in the modes to which these operators are 
applied [Tjj]. The solid lines show the expected theoretical 
dependencies J5), with the vertical scale fit to the data. 



task of process reconstruction is reduced to the known 

HlHliil. This scheme 



problem of state reconstruction 
guarantees that the resulting process is physically consis- 
tent, i.e. completely positive. Furthermore, it permits us 
to incorporate correction for experimental imperfections 
into the reconstruction procedure fl7^ . 

In order to account for the non-deterministic nature of 
the processes being reconstructed, we introduce an ad- 
ditional, fictitious state |0) into the Hilbert space. The 
process can then be treated as deterministic: events in 
which no "click" occurred are interpreted as events in 
which the process has generated state |0) in the output 



11, 22 



The diagonal elements £ J^J 71 of the reconstructed pro- 
cess tensors are shown in Fig. 2J These elements permit 
straightforward interpretation: they give the probability 
that the output of the quantum "black box" contains k 
photons when the m-photon state is present at the in- 
put. According to Eq. (Ql, we expect (fa)^" 1 = mSk, m -i 
and (£a)kk n = (m + l)Sk, m +l- The experimental result 
is consistent with this expectation and explicitly features 
the "raising" and "lowering" properties of and a. For 
input state |m), the output state is similar to |m + 1) for 
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FIG. 4. The "diagonal" values of the process tensor E™™ 
of the photon annihilation (a) and and creation (b) recon- 
structed from the experimental data, with correction for ex- 
perimental imperfections. Each color corresponds to the pho- 
ton number distribution in the output state for the Fock state 
\m) at the input. Insets: worst-case fidelities of the recon- 
structed processes acting within subspaces H' & at (n) of the 
optical Hilbert space spanned by number states |1) . . . |n) and 
0) ... | n — 1), for a and , respectively. 



operator and \m — 1) for operator a . The height of 
the bars in the plots increases linearly with m, which is 
associated with the squared factors \Jm + 1 and y/rn in 
the right-hand sides of Eq. ([1]). 

The consistency of the estimated process tensors with 
those theoretically expected can be quantified using the 
fidelity benchmark. We estimate the worst-case fidelity 



between normalized states a I 



,a< 



and the respec- 



tive outputs fa, at (IV') {^PD °f the reconstructed processes. 
The processes are applied to pure states within subspaces 
H'. - t (n) of the optical Hilbert space spanned by number 
states |1) . . . \n) and |0) . . . \n — 1), respectively, with n 
ranging between 1 and n max . For each n, we employ the 
genetic algorithm to identify the input state producing 



the lowest fidelity in H' & _ t (n) and plot the correspond- 
ing fidelity in the insets of Fig. [4] The fidelities decrease 
somewhat with increasing n because, for high photon- 
number states, the overlap with the probe states is low 
and hence the experimental data do not provide suffi- 
cient information about the effect of the process on these 
states 
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To summarize, we have experimentally implemented 
the operations of optical photon creation and annihila- 
tion and used the technique of coherent-state quantum- 
process tomography to explicitly evaluate, for the first 
time, their process tensors. The reconstructed process 
tensors exhibit "raising" and "lowering" properties of 
these operations. This is the first experiment in which 
complete tomography of trace-non-preserving quantum 
processes has been carried out. 

We thank B. Sanders and C. Simon for helpful dis- 
cussions and A. J. Hendricks, A. Chandra, A. Fedorov 
and A. S. Prasad for help in the lab. We acknowledge 
financial support from NSERC and CIFAR. 



SUPPLEMENTARY INFORMATION 

Theoretical background. Here we explain why the se- 
tups shown in Fig. [T] approximately implement the action 
of the photon annihilation and creation operators. The 
action of a beam splitter on quantum states in the sig- 
nal and trigger input modes described by annihilation 
operators a s and a t can be modeled as evolution under 
Hamiltonian 125] 



H, 



BS 



i\a s a\ — i\a\a t , 



(6) 



where A, assumed real, is related to the amplitude reflec- 
tivity r of the beam splitter as r = sin(£), with £ = Xr/h 
and t being the fictitious interaction time. If |£| <C 1, the 
evolution under Hamiltonian (|6]) can be approximated as 



j-iHr/h 



1 - iHr/h. 



(7) 



The beam splitter input consists of the target state 
in the signal channel and the vacuum state in the trigger 
channel. The evolution operator, acting on this state, 
yields 



-iffBST/S 



(Wjo> t )»|V>.|o) t + c(a 4 |V).)|i> t . (8) 



Conditioning on single-photon detection in the trigger 
channel, we obtain state a s \tp) s in the signal. 

The parametric down-conversion is characterized by 
the simplified Hamiltonian 



PDC 



= i\a\a\ — iXa. 



a i . 



(9) 



In this case, C — Ar/fi has the meaning of the squeezing 
parameter. Applying this Hamiltonian under approxima- 
tion (J7J to the tensor product of the target and vacuum 
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state, we obtain 



G^)jO) t )«|V)JO) t +C(at 



(10) 



which projects onto the signal state a\ \ip) s if a trigger 
photon is detected. 

The above treatment is valid as long as higher-order 
terms in the Taylor decomposition (f[J) can be neglected. 
This is the case when C 71 !^) ^ 1> with nw\ being the 
mean photon number in the target state. 

Experimental imperfections. The models used to ac- 
count for the experimental imperfections are shown in 
Fig. [5] For the process associated with operator a, 
the detrimental effects are the linear losses, non-unitary 
quantum efficiency and the electronic noise of the homo- 
dyne detector |26|. All these effects can be quantified 
and their cumulative contribution modeled by an atten- 
uator with transmission T\ — 0.75 placed after a "black 
box" containing an ideal photon annihilation operator 
[Fig. EJa)]. Accounting for linear losses in maximum- 
likelihood homodyne reconstruction is a well-known tech- 
nique, which consists of modifying the measurement op- 
erator associated with detecting field quadrature values 

[HEllSi. 
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FIG. 5. Model of experimental imperfections for the photon 
annihilation (a) and creation (b) setups. The top and bottom 
schemes in (b) are equivalent to each other. Notation: BS, 
beam splitter; A, attenuator; the quantities in parentheses 
denote the transmission of the corresponding optical element. 
D denotes the operator of phase-space displacement by AX = 
\/2Ti (1 — T^Qin- The amplitude of the input coherent probe 
state is a- m while a is the effective amplitude of the coherent 
state in the mode acted upon by the operators. 

For the photon creation, correcting for the imperfec- 
tions is more complicated, because one needs to take into 
account the mode mismatch between the probe field and 
the signal mode of parametric down-conversion, which is 
determined by the mode of the pump and the optics in 
the idler channel [27|- This mismatch is modeled by a 



Mach-Zehnder interferometer with beam splitting ratio 
Ti and an ideal photon creation operator placed into one 
of its arms [Fig. [SJb), top]. To determine T 2 , we ob- 
serve that, when the input |ai n ) is in the vacuum state, 
the output is expected to be a statistical mixture of the 
single-photon and vacuum states, with the single-photon 
fraction equal to T\T 2 . The output that we observe ex- 
perimentally in this case is indeed, to a high degree of 

with the 
0.75, 



precision, described by such a mixture [28 
single-photon fraction of 0.59. Knowing that T\ 
we conclude that T 2 = 0.59/0.75 = 0.79. 



With this model, if the coherent state entering the in- 
terferometer is |ai n ), the state entering the "black box" 
is | a) = |%/X2ai n ). The latter is the amplitude that we 
refer to in Figs.[2]and[3][b) as well as when using Eq. ([3]) 
to evaluate the process tensor of the photon creation op- 
erator [Fig. Hfb)]. Note that our model is corroborated 
by the photon statistics shown in Fig. EJb). The fit to 
the experimental data yields pr at (a) oc 1 + 0.98a 2 which 
is close to the theoretically expected pr at (a) a 1 + a 2 . 

For the purpose of accounting for the imperfections 
in the process reconstruction, it is convenient to re- 
formulate their model in terms of the scheme shown 
in the bottom panel of Fig. [5jb). The output of the 
"black box" undergoes a linear loss channel with trans- 
missivity T\Ti followed by phase-space displacement by 
AA = v / 27Y(l — T2)u\ a . These transformations must be 
compensated for in the process reconstruction. To that 
end, prior to launching the iterative algorithm with cor- 
rection for a linear loss, we apply reverse displacement to 
the measured quadrature data. Specifically, we subtract 
AA cos 8 from each experimentally measured sample of 
quadrature observable Xg — X cos 9 + P sin 8, where 9 is 
the local oscillator phase. The modified values are then 
used as input for the iterative algorithm. 

Calibrating coherent state amplitudes. The probe co- 
herent states are generated by a series of three vari- 
able attenuators, each made up of a half-wave plate 
and a polarizer, followed by neutral-density attenuators. 
The amplitude of the coherent state is then given by 
a- m = A sin 29i sin 28 2 sin 28%, where 8's are the waveplate 
orientation angles and A is a common factor. Determin- 
ing A directly from power measurements proved unprac- 
tical. Instead, we calibrated this factor by measuring the 
amplitudes of the coherent states |a ou t) observed with 
the homodyne detector in the absence of trigger events 
in the photon annihilation setup [Fig. Ufa)]. This ap- 
proach is more reliable because an absolute gauge for 
the amplitude scale is provided by the quadrature vari- 
ance of the vacuum state. Once a ut has been evaluated 
(Fig. [6]), the input coherent state amplitude is determined 
by taking into account detection losses, i.e. according to 

«out = VTl^in- 
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amplitude according to waveplate angles 
sin(29 1 )sin(2e 2 )sin(2e 3 ) 

FIG. 6. Calibration of the coherent state amplitude by ho- 
modyne measurement. 
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